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ABSTRACT 

We study gaseous outflows from disk galaxies driven by radiation pressure on dust grains. We include 
the effect of bulge and dark matter halo and show that the existence of such an outflow implies a 
maximum value of disk mass-to-light ratio. We show that the terminal wind speed is proportional 
to the disk rotation speed in the limit of a cold gaseous outflow, and that in general there is a 
contribution from the gas sound speed. Using the mean opacity of dust grains and the evolution of 
the luminosity of a simple stellar population, we then show that the ratio of the wind terminal speed 
(voo) to the galaxy rotation speed (y c ) ranges between 2-3 for a period of ~ 10 Myr after a burst of 
star formation, after which it rapidly decays. This result is independent of any free parameter and 
depends only on the luminosity of the stellar population and on the relation between disk and dark 
matter halo parameters. We briefly discuss the possible implications of our results. 
Subject headings: galaxies: starburst — galaxies: evolution — intergalactic medium 



1. INTRODUCTION 

Galactic outflows are believed to be an important in- 
gredient in the evolution of galaxies. It is thought to 
provide a feedback mechanism for the regulation of star 
formation in galaxies, and the subsequent chemical evo- 
lution in them. Also, galactic scale winds are thought 
to enrich the intergalactic medium (IGM) with metals. 
Observations also show that galactic scale outflows are 
common in starburst galaxies in both local and high red- 
shift universe (Veilleux et al. 2005). 

In the standard scenario, the interstellar medium 
(ISM) of the starburst galaxy is heated by supernovae 
(SN) , the thermal pressure of the hot gas drives the out- 
flowing gas. Larson (1974), Saito (1979), Dekel & Silk 
(1986) argued that supernovae in a star forming galaxy 
would drive an outflow in excess of the escape speed from 
the galaxy. This scenario, however, has met with prob- 
lems from new observations. Recent observations of cold 
clouds embedded in the hot gas show that the maximum 
speed of these clouds is correlated with the star forma- 
tion rate (Martin 2005), a correlation that is more easily 
explained by outflows driven by radiation pressure than 
thermal pressure. Observations of Lyman break galaxies 
(LBG) at z ~ 3 have also found a correlation between 
the outflow speed and the star formation rate, as well as 
with the reddening due to dust. Also, SN explosions in 
disk galaxies may only produce a thickening of the disk 
gas because of the incoherent nature of these explosions 
(Fragile, Murray, Lin 2004). Moreover, cold clouds em- 
bedded in the hot gas in the outflow are not likely to 
survive longer than a Myr, because of various instabili- 
ties and/or evaporation by thermal conduction. 

Murray et al. (2005; hereafter referred to as MQT05) 
showed that radiation pressure can be comparable to the 
ram pressure in hot outflows, and considered radiatively 
driven shells of gas and dust. Martin (2005) found that 
this scenario provides a better framework in order to un- 
derstand the observed correlations. Nath & Silk (2008) 
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discussed a hybrid model of outflows with radiation and 
thermal pressure. 

In this paper, we study dust driven gaseous winds from 
luminous disk galaxies. We estimate the terminal speed 
of the combined gas and dust flow, taking into account 
the gravity of the dark matter halo, and the disk mass-to- 
light ratio that is expected from stellar population syn- 
thesis models. 

2. DUST-DRIVEN WIND FROM A KEPLERIAN DISK 

Recently Zhang & Thompson (2010) considered a disk 
of radius d constant surface density (£) and surface 
brightness (/). Assuming a cylindrical geometry, the 
force of gravitation f g (z), and that due to radiation 
f r (z), along the pole are given by, 



fg = 27TGS 
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(r 2 + z 2 ) 3 / 2 
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where n is the average opacity of dust grains to absorp- 
tion and scattering of photons. The ratio of these forces, 
the Eddington ratio, increases with the height z, begin- 
ning with a value of Tq = 2 cGT, a ^ * nc disk centre at 
z = 0. From energy conservation they found the velocity 
at height z, to be, 
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This implies a terminal velocity of v c 
\/47rG£<i(7rro/2 — 1), with a lower limit on Tq 



(2) 
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for driving an outflow. 

Assuming a Keplerian disk, with average rotation 
speed of (v C: kep) ~ \ V TrGTid, the terminal velocity can 

be written as = ^{\J ^f- — l)(v c ,kep), which yields 

a value Voo — 3(v c ,kep), for To = 1. But these results 
will change with the proper inclusion of bulge and dark 
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matter halo. 



3. GASEOUS OUTFLOWS 



In this paper, we discuss steady, rotating wind in which 
dust is propelled outwards by radiation pressure and 
drags the gas with it. We ignore magnetic forces and 
treat the wind as a single-phase fluid, which is marginally 
optically thick (MQT05). 

Although our results are valid for all temperatures, in 
practice the model allows for only cold flows (as observed 
in Nal and Mgll lines (e.g. Martin 2005)) for the fol- 
lowing reasons. The sputtering radius for grains over 
a wind timescale of 10 Myr (which we derive below) is 
a ~ 50 x oT7S A for T - 10 5 K (Tielens et al. 1994), 

which implies that only grains smaller than 5 A are de- 
stroyed, without changing the opacity. But the radiative 
cooling time is small (< 1 Myr), so for adiabatic ap- 
proximation to be valid, the adiabatic-cooling time scale 
should be shorter than the radiative-cooling timescale, 
which limits the temperature < 10 4 K. Also the timescale 
for dust-heating is ~ 10 Myr and can be neglected. 

Secondly, we improve upon the previous estimates of 
the required value of Tq needed to drive a wind, by tak- 
ing into account the gravity of bulge and halo. Clearly 
the relation between Voo of the wind and v c in which 
the rotation speed is estimated from a keplerian disk is 
not relevant for disk galaxies, because the rotation speed 
must be calculated from the dark matter halo. 



3.2. Wind terminal speed 

To determine the terminal speed of the wind, we use 
the fact that the Bernoulli function is preserved along 
a streamline, assuming that a streamline extends from 
the base to infinity. In an isothermal wind the terminal 
speed tends to infinity as the wind maintains constant 
sound speed. It is however more reasonable to assume a 
polytropic equation of state. 

One can write the Bernoulli equation for a polytropic 

2 

gas (with adiabatic index 7) along a streamline: ^ + 

2 

+ <f> = E, where c s is the sound speed, <f> is the 
potential and E is a constant. Equating the values at 
the base and infinity, we have 
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where Vf, is the wind speed at the base (z = 0), c S; h is the 
sound speed at the base, and c Sj00 is the sound speed at 
infinity, which is negligible. The potential is, 
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3.1. Disk, bulge and halo parameters 

We consider a disk with constant surface density (E) 
and surface brightness (I), of radius d, and which is em- 
bedded in a bulge and a halo. We assume a spherical 
mass distribution in the bulge and the halo. For the 
bulge, we assume a total mass of M& inside a radius 
n < d. 

For the halo, we consider a Navarro-Frenk- White 
(NFW) profile, with total mass Mh characterized by a 
concentration parameter c = r V i r /r s (Navarro, Frenk, 
White 1997). We fix the total halo mass for a given 
disk mass (M d = 7rd 2 E), by the ratio M h /M d ~ 1/0.05, 
as determined by Mo, Mao & White (1998) (referred 
to as MMW98 hereafter). We use the prescription of 
MMW98 for the disk exponential scale-length R d ~ 

(l/v / 2)Ai?2oo/<T 1 ^ 2 /fl, using their eqns 23, 26, 32, and 
use d = Rd/V2, since the total masses in the case 
of uniform density and exponential disk are given by 
M = 7r<i 2 E = 27rEoi? 2 >. The rotation speed implied by 
the NFW profile peaks at a radius r ~ 2r s , given by, 



ln(3) - 2/3 



■200 



2 ln(l + c) - c/(l + c) 



(3) 



where r s is the scale radius of NFW profile and i>2oo is the 
rotation speed at the virial radius. We choose this value 
of the maximum rotation speed to represent the v c of the 
disk galaxy, since Figure 2 of MMW98 shows that the 
value of v c from the flat part of the total rotation curve 
does not differ much from the peak of the rotation curve 
from halo only. We also use r b /d ~ 0.1, and M b /Md ~ 
0.5, consistent with observed range of luminosity ratio 
between bulge and disk (Binney & Mcrrificld 1998). 



Here the first term denotes a pseudo-potential due to 
radiation pressure. The second term refers to the gravi- 
tational potential of the disk, and here we have assumed 
for analytical simplicity that gas stays near the pole, so 
these forces are given by eqn([T|). The third term is the 
potential due to the outward centrifugal force of the ro- 
tating gas in the wind, and the last two terms denote the 
effect of the bulge and halo gravity. We also assume that 
the bulge exerts a negligible radiation pressure, since the 
dominant bulge stellar population is old and red and the 
mean opacity k of dust grains in these wavelengths is 
smaller than in blue band. 

The centrifugal force at the base of the rotating wind 
should be equal to gravitational force due to the mass 
inside a, the radial distance at the base. We assume that 
a > 7'b, the bulge radius, while being much smaller than 
d, the disk radius. Hence the centrifugal force is written 



Hase GM b 
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The first term on RHS is due to the bulge, and the second 
term, due to the dark matter halo. For a <C r s , the 
second term can be approximated as G 2 ^{ s which is much 

smaller than the bulge term ^2, because the gravity 
of baryons dominates at small radii. Neglecting the halo 
term, the specific angular momentum is given by, 



y/GM b a . 



(7) 



Finally, the values of the potential at the base and infinity 
are as follows, 
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Fig. 1. — The ratio between Voo and v c is shown as a function 
of To for galaxies with two different concentration parameter c, for 
cold gaseous outflow. 
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We note that the radiation pseudo-potential is zero at 
the base but ~ —2-kGYjoHJ-^-) at infinity along the pole. 
Since the dark matter halo is truncated at r200; its po- 
tential vanishes at infinity. 

It is reasonable to assume that the wind speed at 
the base is comparable to the sound speed in the disk 
(vb ~ c s .b). Putting it all in cqn(0|, we have the follow- 
ing expression for the terminal speed, 
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Wc emphasize that Tq refers only to disk parameters and 
To = 1 does not signify Eddington luminosity for radia- 
tion pressure on dust. 

Figure 1 shows the dependence of Voo/v c on To = 
kI /2-kcGTi, and on the concentration parameter c for 
c St b = 0, and for a = r b . The sound speed makes only a 
small difference to the terminal speed, but it helps estab- 
lishing a quasi-steady regime in the outflow redistribut- 
ing the pressure over the entire flow. The curves show 
that there is a minimum value of Tq ~ 2 for the wind to 
reach infinity. In general the wind speeds are lower for 
higher c, which is the case of low-mass and more com- 
pact galaxies, but the variation is small. The important 
finding here is that the wind requires To > 1. 

We next use a relation between c and galactic mass 
given by Maccio et al. (2007), and calculate the wind 
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v o =200 km/s,z = 
v c =100 km/s,z = 
v c =200 km/s,z=7 



Fig. 2. — The ratio between Vaa and v c is shown as a function of 
To for galaxies with different v c at two different redshifts. 

speed for galaxies with different masses, or rotation 
speeds. One can extend the calculations to higher red- 
shift by using the scaling c oc (l + z) -1 . We use 
Q = 0.27, h = 0.7, and calculate the NFW profiles for 
the A-cosmology. Figure 2 shows the results for Voo/v c 
as a function of Tq for different v c , at z = and at z = 7. 
The curves show that the v^jvc depends mostly on To 
and varies weakly with rotation speed and rcdshift. 

We next ask the question what determines the value of 
Tq which is crucial for the calculation of v^. 



3.3. Evolution of wind 



with time 



We recall that the value of To depends on disk param- 
eters (£,-0 and dust grain properties (through k). Con- 
sider first the value of £// that is essentially the disk 
mass-to-light ratio. One can compare the observed val- 
ues of mass-to-light ratio with the minimum requirement 
as derived above. Li and Draine (2001) give the mean 
opacity for gas mixed with dust as ~ 128 cm 2 g" 1 in U 
and ~ 93 cm 2 g _1 in the B band. Here wc shall con- 
sider the B band value as a conservative estimate, which 
ensures that star formation is not obscured. In case of 
obscured star formation, shell(s) of gas and dust would 
be accelerated, which are likely to undergo fragmentation 
and produce clouds due to instabilities, facilitating the 
escape of ionizing photons (Razoumov & Sommer-Larsen 
2010; Wise & Cen 2009). 

Using B-band values we find a maximum value of disk 
mass-to-light ratio required for the outflow to occur, 
given by M/L ~ k/ (2cGr ). For c ~ 10, this is given by 

,$ M/ u Ie \ = Tb < 0.036. This upper limit on disk T B 

(Lb/Lb, o) " — ^ L " 

is much lower than observed in present day disk galaxies. 
The Milky Way disk has a local value of T B ~ 1.2 ± 0.2 
(Flynn et al. 2006), and the typical value for disk galaxies 
is 1.5 ±0.4 (Fukugita, Hogan, Peebles 1998). According 
to Flynn et al. (2006), a third of the total disk mass 
comes from gas and the rest from stars. Starburst galax- 
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ies can have a much lower value of Tb- The estimated 
M/L bol for NGC 7714 is - 0.02 (Bernlohr 1993a), for re- 
gions in NGC 520 it is - 0.003 (Bernlohr 1993), and for 
a young super cluster in M82, Smith & Gallagher (2001) 
has estimated the Yy ~ 0.02. These low values of T are 
believed to arise from a top heavy IMF and young age of 
the stellar population (e.g., Kotilainen et al. 2001; Smith 
& Gallagher 2001). 

For an instantaneous burst these models predict an ini- 
tial period of roughly constant luminosity for t < 3 Myr, 
and a decrease in the luminosity afterwards (Bruzual & 
Chariot 2005; Buzzoni 2005; Vasquez & Leitherer 2005). 
For a Salpcter IMF and a stellar mass-luminosity relation 
of the type L <x M° , the late-time decay of the luminos- 
ity is given by L cx t-^- 1 - 35 )/^- 1 )^ ^°' 9 , for /3 = 3.5). 
The initial period of rather constant luminosity stems 
from the fact that while low-mass stars are yet to col- 
lapse, the massive stars evolve quickly, and the duration 
of this period corresponds to the main-sequence life-time 
of the most massive stars. 

Using these models we can determine the T for a stellar 
population, multiplying by a factor ~ 3/2 to account for 
an additional gas mass, and determine the expected disk 
Tb ratio as a function of time after an instantaneous 
starburst. Using the results in Figure 9 of Vasquez & 
Leitherer (2005), which uses a Salpeter IMF between 0.1 
and 100 M , and using Mb,q = 5.45, one can calculate 
T b as a function of time. Strength of instantaneous star- 
burst is characterized by total mass converted into stars 
initially, which in this case is 10 6 M Q . 

Using the mean dust opacity for B band, we calculated 
the time evolution of To, and then determined the wind 
terminal speed using cqn®. Although the use of eqn® 
assumes a constant Tq along a streamline, the distance 
travelled by the wind over the time-scale of change of Tq 
(~ lOMyr) is large (I > 6.3kpc(u/600km/s)). In other 
words, the time-scale for the wind to reach a considerable 
height above the disk is comparable to the time-scale of 
evolution in To, and therefore we can use our formalism 
to estimate the terminal speed with the evolution in Tq. 

Figure 3 shows the evolution of Wqo/^c with time for 
z = (solid line) and z = 7 (dotted line). The curves 
show that the wind speed decreases rapidly after ~ 10 
Myr, and that Voo < 3v c , its value being smaller for com- 
pact galaxies. This result can be compared with the ob- 
served range of maximum wind speed. Martin (2005) 
found that the maximum speed of clouds embedded in 
outflowing gas ranges between 2-3 v c , and Rupkc ct al. 
(2005) found a range of [0.67-3] v c . We show this range 
with two dashed lines in Fig 3. Figure 3 also shows that 
the wind speed is somewhat smaller at high redshift. The 
reason is that galactic mass for a given v c is smaller at 
high redshift, but c oc M~ 0,2 (1 + z) , and the mass 
effect outweighs the redshift effect. However, the varia- 
tion of Voo/vc with redshift is expected to be very small 
in this model, much smaller than those caused by other 
parameters, such as the IMF. 

4. DISCUSSIONS 

Figure 3 shows that dust-driven winds are likely to 
have a terminal speed ~ 2-3i> c , for a combination of rea- 
sons that involve stellar physics and the relation between 
disk and halo parameters. It is interesting that this result 
coincides with observations, since there is no free param- 
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Fig. 3. — The evolution of terminal wind speed in terms of rota- 
tion speed with time is plotted for v c = 100 km/s at z = (solid 
line), and at z = 7 (dotted line). Horizontal dashed lines show the 
observed range at z = (Rupke et al. 2005, Martin 2005). 

eter in our calculation. The strength of our approach lies 
in the fact that the terminal speed calculated using the 
Bernoulli function is independent of the streamline used 
by the gas, as long as streamlines do extend to infinity, 
which is our basic assumption. Below we discuss a few 
implications. 

In the scenario of energy driven winds, the IGM is be- 
lieved to be enriched by winds from dwarf galaxies, since 
they were numerous in the early universe (Silk, Wyse & 
Shields 1987; Nath & Trentham 1997; Ferrara, Pettini & 
Shchekinov 2000; Cen & Bryan 2001; Madau, Ferrara, 
Rees 2001; Aguirre ct al. 2001). However, in the case of 
dust-driven winds, the importance of low-mass galaxies 
in IGM enrichment diminishes because oc v c . Our cal- 
culations here show that the wind speed depends strongly 
on the time elapsed after a starburst, or more generally 
on the star formation history and parameters. It is be- 
lieved that the IMF is weighted towards massive stars at 
high redshift (e.g., Schneider & Omukai 2010), in which 
case the wind speed likely increases with redshift. In this 
case the contribution of dwarf galaxies at high redshift 
may still be important. 

Recent simulations for IGM enrichment including mo- 
mentum driven winds have used the ansatz Voo ~ 3v c 
(Oppcnheimcr & Dave 2006), using the estimates of 
MQT05 for wind from a spherical galaxy radiating close 
to the dust Eddington luminosity. It has been suggested 
that the momentum driven winds drive a feedback loop 
that makes v^, ~ 3v c (Martin 2005; MQT05). Our cal- 
culations show that such dust-driven winds are possible 
only within a period of ~ 10 Myr after a starburst. Since 
the wind speed depends strongly on the value of disk T, 
which depends on the IMF (being smaller for a top-heavy 
IMF), the wind speed is expected to be larger for a IMF 
skewed towards massive stars. The implications needs 
to be studied with numerical simulations using modified 
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prescriptions for the wind speed. It is possible that this 
puts constraints on the metallicity of the wind gas. In 
addition, as the dust opacity is proportional to the metal- 
licity, k oc Z, the existence of a threshold Tq needed for 
the wind to reach infinity suggests that only metal-rich 
galaxies can enrich IGM. 

We note that our calculation does not determine 
streamlines, without which we cannot calculate the den- 
sity structure in the wind, and therefore cannot derive 
the mass loss rate. In the case of single scattering, the 
maximum mass loss rate is M ~ L/(vooc) (MQT05). 
This recovers the scaling that M oc vj 1 adopted by nu- 
merical simulations (Oppcnheimer & Dave 2006). Taking 
our results into account, the mass loss rate for a given 
luminosity is expected to be lower at high redshift, be- 



cause of the possible rise of with z arising from IMF 
evolution. 

In summary, we have derived a terminal speed for dust- 
driven outflows from disk galaxies, and have shown that 
Woo ~ (2-3)v c , which is determined by the minimum 
value of disk ~ 10~ 2 arising from the luminosity 
of a stellar population, and the relation between the disk 
and dark matter halo that fixes the terminal speed for a 
given value of T. We have shown that dust-driven winds 
from disk galaxies are excited within a time-scale of ~ 10 
Myr of a starburst, after which the radiation pressure on 
dust is unable to drive outflows. 

We thank Mitchel C. Begelman, J. Silk and S. Srid- 
har for valuable discussions. We thank the referee CD. 
Matzncr for useful comments. 
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